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Adaptive models ol opinion formation among humans can display a fragmentation transition, 
where a social network breaks into disconnected components. Here, we investigate this transition in 
a class of models with arbitrary number of opinions. In contrast to previous work we do not assume 
that opinions are equidistant or arranged on a one-dimensional conceptual axis. Our investigation 
reveals detailed analytical results on fragmentations in a three-opinion model, which are confirmed 
by agent-based simulations. Furthermore, we show that in certain models the number of opinions 
can be reduced without affecting the fragmentation points. 



I. INTRODUCTION 

In many different fields networks have been used to 
describe and analyze complex systems consisting of in- 
teracting subunits. The applications of networks range 
from biological systems to technical devices and social 
communities PHI] ■ Accordingly, the building blocks of a 
network, the network nodes, can correspond to different 
entities, such as genes, neurons, computers, websites or 
individuals. The interactions among them, the links be- 
tween the nodes, represent e.g. chemical reactions, phys- 
ical connections, or social bonds. In the applications the 
temporal evolution of a network is often governed by two 
different types of dynamics: the dynamics on the net- 
work, describing the evolution of the internal degrees of 
freedom, and the dynamics of the network, capturing the 
evolution of the network topology. 

Adaptive networks are characterized by an interplay 
of the dynamics on the network and the dynamics of the 
network, where neither of both types of dynamics can be 
neglected |5] . It has been shown that this interplay 
gives rise to the emergence of complex topologies and 
dynamics [7] , spontaneous appearance of different classes 
of nodes from an initially homogeneous population [5J [S] , 
and robust self-organization to critical states associated 
with phase transitions [10l HI] • The self-organization of 
adaptive networks is believed to be of importance in the 
evolution of cooperation, opinion formation processes, 
epidemic dynamics, neural networks, and gene regulation 

In the adaptive networks literature, opinion dynam- 
ics has currently attracted particular attention [13H18) . 
Typically, in these models a society is described as a net- 
work, where nodes correspond to individuals and links 
to social relationships. The internal state of the indi- 
viduals indicate their position regarding some issue, such 
as political opinion, religious affiliation, or musical taste. 
The individuals change their states by adopting opinions 
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from their topological neighbors. The network topology, 
i.e. the specific pattern of nodes and links, changes as 
individuals break up relationships with dissenting neigh- 
bors and/or establish new relationships to those holding 
similar opinions. 

The simplest adaptive-network model for opinion for- 
mation is the adaptive voter model [13j IT9H2T] . In this 
model and in many of its variants [TJ1 the relative 

rate of change of the topology compared to the change 
of node states is controlled by a single parameter, p, the 
rewiring rate. Depending on this parameter, the net- 
work either reaches a global homogeneous state, where 
all nodes hold the same opinion or a fragmented state, 
consisting of two disconnected components which are in- 
ternally homogeneous. The transition separating these 
two regimes is called a fragmentation transition. 

While most of the voter-like models only consider a bi- 
nary choice of opinions, many real world situations offer 
a larger number of choices. In the physics literature some 
models for opinion formation, which consider arbitrary- 
many opinions have been studied [TJ1 H3] . In these mod- 
els all opinions are "equidistant" in the sense that all in- 
teractions between any given pair of (different) opinions 
follow the same dynamical rules. Models recognizing that 
the outcome of interactions may depend on a measure of 
similarity (or distance) between opinions are often con- 
sidering an uncountable set of opinions and are therefore 
hard to treat analytically [T5J [55] . For instance in [55] 
opinions are placed on a 1-dimensional axis, resembling 
e.g. the political spectrum. 

Here, we consider a natural extension of the original 
adaptive voter model, where we allow for an arbitrary 
countable set of opinions. In the proposed model the 
rewiring rate that governs the interaction of conflicting 
agents is assumed to depend on the specific pairing of 
opinions held by the agents. The model can thus ac- 
count for heterogeneous "distances" between opinions. A 
large distance, characterized by a high rewiring rate, in- 
dicates a controversial pairing, whereas a small distance, 
and correspondingly low rewiring rate, indicates that the 
respective opinions are almost in agreement. 

The proposed model is described in detail in Sec.|TI] We 



2 



then calculate the fragmentation diagram of a three-state 



model in Sec. Ill and derive the corresponding fragmen- 
tation thresholds for systems with an arbitrary number 
of states in Sec. |IV| Finally, in Sec. |Vj we show that in 
certain systems a reduction is possible such that the dy- 
namics can be captured by a system with a lower number 
of opinions. 



II. MULTI-STATE VOTER MODEL 




We consider a network of N nodes, corresponding to 
individuals, and L links, corresponding to social con- 
tacts. Each node a holds a state s a , indicating the 
opinion held by the corresponding individual. The net- 
work is initialized as a random graph with mean degree 
(k) = 2L/N. The initial node states are drawn ran- 
domly and with equal probability from the set of all states 
r = {gi, c/2, 53, ■ ■ • j 9g}, where the total number of states 
is |r| = G < N. 

The system is then updated as follows: In each update 
step a random link (ct,/3) is chosen. If s a — sp the link 
is said to be inert and nothing happens. If s a =/= sp, 
the link is said to be active and an update occurs on 
the link. A given update is either a rewiring event or an 
opinion adoption event, decided randomly depending on 
the similarity of the respective opinions. For individuals 
a, P with opinions s a = gi and sp — gj, the update 
is an rewiring event with probability and an opinion 
adoption event otherwise (probability 1 — Pij). In the 
following the parameters ptj are called rewiring rates. 

In a rewiring event the focal link (a, /3) is severed, and 
a new link is created either from a to a randomly chosen 
node 7 with s 7 = s a , or from /3 to a randomly chosen 
node 7 with s 7 = sp. The choice between the two out- 
comes is made randomly with equal probability. In an 
opinion update, either node a changes its state to sp 
or node j3 changes its state to s a , where the choice be- 
tween both outcomes is again made randomly with equal 
probability. In the following we assume symmetric in- 
teractions, which implies Pij = pji such that the specific 
model is characterized by a set {p^} of G(G — l)/2 pa- 
rameters. 

The model proposed above preserves the symmetry 
of the direct interaction of two opinions postulated in 
the adaptive voter model, i.e. in direct comparison no 
opinion is stronger then the other. However, it breaks 
the symmetry between different pairings of opinions such 
that rewiring is more likely in certain pairings than in 
others. 

From the adaptive voter model [131 119] it is known 
that the fragmentation transition separates a so-called 
active regime from a fragmented regime. In the active 
regime a finite density of active links persists in the long- 
term dynamics such that there is ongoing dynamics until 
fluctuations drive the system eventually to an absorbing 
consensus state. In the fragmented regime, disconnected 
components emerge, which are internally in consensus. 



FIG. 1. State-network of a three-state voter model. The 
nodes in this network are the three opinions A, B and C. 
The rewiring rates pi,P2 or p$, encode the degree of contro- 
versy between agents holding the respective pairings of opin- 
ions AB, BC or CA. The ellipse and dashed lines illustrate 
the example described in the text: partial fragmentation with 
respect to A, i.e. the fragmentation of the corresponding net- 
works of agents in a component where all individuals hold 
opinion A and a second component where the individuals hold 
opinions B and C. The ongoing dynamics in the latter com- 
ponent may later lead to the disappearance of either opinion 
B otC. 



Because the G-state model contains several different 
types of active links (corresponding to all possible pair- 
ings in r), regimes can occur where active links of a cer- 
tain type vanish while others prevail. This can lead to 
configurations where a certain subset of the states only 
appears in one component of the network in which no 
state not belonging to this subset is present. In the fol- 
lowing we call this situation a partially fragmented state. 
In contrast to the fully fragmented state where every 
component is internally in consensus, the dynamics in 
the partially fragmented state can continue in some com- 
ponents while others may be frozen in internal consen- 
sus. The partial fragmentation cannot be undone, so 
that achieving global consensus is impossible after par- 
tial fragmentation has occured. However, the ongoing 
dynamics in the active components will eventually lead 
to internal consensus in every component. The absorbing 
state which is ultimately reached after a partial fragmen- 
tation therefore consists of 1 < 7 < G major components, 
holding the 7 surviving opinions, respectively. For 7 = G 
we recover full fragmentation and the case where 7 = 1 
we denote as the fully active regime where all types of 
active links prevail. Only in the latter case (due to finite- 
size effects) global consensus can be reached eventually. 

III. FRAGMENTATION TRANSITIONS IN A 
THREE-STATE VOTER MODEL 

We start our exploration of the proposed multi-state 
voter model by considering the case G = 3, which is 
the simplest case which is not trivial (G = 1) or exten- 
sively studied (G = 2) jTBl EH |25]. Let us consider the 
set of opinions T = {A, B,C}, giving rise to three dif- 
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ferent rewiring rates {pab,Pac,Pbc} which we denote 
as {pi,P2,P3} according to the state-network depicted in 

Fig.g 

Let us emphasize again the difference between the 
state-network and the network of individuals. The 
state-network is a complete, weighted graph with G 
nodes which represents the relationships between differ- 
ent states, such that states connected by small rewiring 
rates are similar to each other, whereas states which are 
connected by large rewiring rates differ significantly from 
each other. The network of individuals, in contrast, is an 
unweighted graph with N nodes and mean-degree (fc), 
which represents the interactions between individuals. 

In principle, the three-state system can reach five dif- 
ferent final states: fully active (leading eventually to 
global consensus), full fragmentation, and partial frag- 
mentation with respect to A, B, or C. Here, partial 
fragmentation with respect to a certain state refers to a 
situation where a component of nodes in that particu- 
lar state fragments from an active component (a mixed 
component of nodes in the other states). 

For the calculation of fragmentation thresholds we fol- 
low the approach given in [16] . We determine the evolu- 
tion equations for the number of active motifs (network 
motifs containing active links) starting from a situation 
close to the fragmentation threshold. For simplicity, the 
present article uses only the simpler of two different ac- 
tive motif bases proposed in [TB]. We emphasize that 
all calculations below can also be carried out using more 
elaborate motif bases, but at the price of having to deal 
with considerably larger matrices. 

Following To] , we define g-fans as a bundle of q active 
links of one type, say AB-\mks, connected to a single 
A- or i3-node. We do not account for the number of 
inert links connecting to this focal node. For the sake 
of simplicity we also do not consider mixed active motifs 
containing all three states. We confirmed that effects of 
mixed motifs can be suitably captured by the procedure 
described below. 

We start by calculating the condition for partial frag- 
mentation with respect to A (see Fig. [TJ) . The A-cluster 
fragments from the rest of the network when all AB- 
motifs and all AC-motifs vanish. Because in general 
Pi 7^ P2 , we have to treat AB- and AC-motifs separately. 
We start by considering a network with two almost dis- 
connected clusters, one of which is composed purely of 
A-nodes and the other of B- and C-nodes and then ask 
whether the fragmented state is stable, such that frag- 
mentation is reached, or unstable, such that the system 
avoids fragmentation. 

In the almost fragmented state the expected effect of 
network updates on the active motifs is captured by a 
procedure proposed in [Tjj]. For the case of k — 3 we 
obtain the transitions rules shown in Fig. [2] New ac- 
tive motifs are created when an opinion update occurs. 
We approximate the degree of the focal node k by the 
networks mean degree (fc). Because of the clusters be- 
ing almost-separated the newly formed active motif is a 



FIG. 2. Transitions of AB- and ^C-fans for a degree-regular 
network with k = 3 and equiprobable states for the scenario 
of partial fragmentation with respect to A. Black, white and 
grey nodes correspond to agents holding opinion A, B, and C, 
respectively. The active link densities p\ and pi are assumed 
to vanish close to the partial fragmentation point, whereas pz, 
the densitiy of BC-links can be finite. 



k — 1-fan [16] . This fan can subsequently lose active links 
due to subsequent opinion updates and rewiring events. 
We account for a finite density of active i?C-links, P3, in 
the active component by creating an AC-fan (AS-fan) 
instead of an AB-idM (AC-fan) with probability /03 when 
a new fan is created by a S-node (C-node) adopting opin- 
ion A. 

If we start with equal distribution of states, the rela- 
tion p 3 = [BC]/{k[B]) = [BC]/(k[C]) holds, where [B] 
and [C] denote the numbers of B-nodes and C-nodes, re- 
spectively. Note that P3 differs from the global BC-link 
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fragmentation 



density pf ] = [BC]/L. 

The set of transitions for k = 3 (see Fig. |2| defines a 
dynamical system, describing the time evolution of the 
densities of active motifs close to the partial fragmen- 
tation with respect to A. The stability of the partially 
fragmented state in this system is then governed by the 
block-structured Jacobian, 



J(Pl,P2,P3) 

where 



D pi- X pi(P3) X Pl (p 3 ) 

X P 2 (P3) D p 2 - X P2(P3) 



-1 1(1 -p) 4(1 -p) 



D p = 1 |(l-p)-l 



I(l-P)-Iy 



, (1) 



(2) 



and 



'0 i(l -p)p ON 
X p (p)= ( |(1 -p)p 
,0 0/ 



(3) 



The diagonal blocks in the Jacobian given in Eq. ([!]) 
can be interpreted as "self-interaction" terms, captur- 
ing contributions from the same motif-type, and the off- 
diagonal terms as "exchange" terms, capturing contribu- 
tions from different motif-types. The structure of this 
Jacobian remains unchanged for any partial fragmenta- 
tion of a three-state system, while the matrices ^ and 
([3]) change when the motif set is altered. In particular, 
the dimension of these matrices increases with increasing 
mean degree and/or number of motifs considered. 

In a dynamical system a steady state is stable if all 
eigenvalues of the corresponding Jacobian have negative 
real parts [57] • For the present system this means that 
the fragmented state is stable if all eigenvalues of the Ja- 
cobian J are negative and the fragmentation transition 
occurs as at least one of the eigenvalues acquires a pos- 
itive real part. Therefore, demanding A (J) = 0, where 
A(J) is the leading eigenvalue of J, yields a condition for 
the fragmentation transition, depending on the three pa- 
rameters pi,p2, and p 3 . The phase diagram in Fig. [3] is a 
projection of this fragmentation condition on the p\-p 2 - 
plane for the extreme values of p 3 . 

Let us first consider the case where p 3 = 0, which is 
encountered if p 3 exceeds p c , the fragmentation threshold 
of the adaptive two-state voter model. In this case, X 
becomes zero and the set of eigenvalues of the Jacobian 
J is the conjunction of the eigenvalues of the matrices 
D pi and D P2 . Thus, A(J) is negative iff A(D P1 ) and 
A(D P2 ) are negative. Indeed matrices D P1 and D P2 are 
the Jacobians of the two uncoupled two-state systems A— 
B and A — C. Thus fragmentation of A requires that the 
two-state fragmentation condition is met separately for 
the AB and AC subsystems. In other words, if the links 
between B and C nodes vanish (p 3 = 0), fragmentation 
occurs when both pi > p c and p 2 > p c (see Fig. [3] region 
III). 




no fragmentation/ 
consensus 




Pi 



FIG. 3. Fragmentation diagram for a three-state system. 
Fragmentation with respect to opinion A occurs always in 
the black region and never in the grey region. In the dashed 
region, fragmentation can occur depending on the value of 
p3, the density of BC-links in the SC-component. If pi > P3 
and P2 > P3 then this fragmentation diagram characterizes 
the final state of the whole system. In this case region I coin- 
cides with the global consensus regime and full fragmentation 
is reached iff ps > p c , where p c is the fragmentation threshold 
for the adaptive two-state voter model. 



For studying the case p 3 > we first note that every 
matrix- valued row of J sums to D p . , where i = 1,2. Fol- 
lowing [21], as will be discussed below, A(J) is bounded 
by A(D P1 ) and A(D P2 ). Therefore, fragmentation with 
respect to A is guaranteed when pi > p c and p 2 > p c > 
but can already occur when only either p\ > p c or p 2 > p c 
is satisfied (see Fig. [3] region II) . 

The maximum extension of region II is observed when 
p 3 = 0, the corresponding maximal value of p 3 , p m ax, 
can be determined to good approximation by a moment 
closure approach (see Appendix), yielding 



(4) 



2k 



Solving the condition X(3(pi,p2, Pmax)) = numerically 
yields the curve separating regions I and II in Fig. [3] 
Moreover, from the diagram in Fig. [3] it is clear that 
this curve implies the existence of a minimal rewiring 
rate p min , such that for p\ < p min or p 2 < p m \n partial 
fragmentation with respect to A becomes impossible. 

Let us emphasize, that calculations of fragmentation 
thresholds for partial fragmentations build on the esti- 
mation of the active link densities from given rewiring 
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FIG. 4. Numerical phase diagram for the three-state model. 
Color-coded is the density of links, connecting the A and the 
BC-cluster (pi + P2) over the rewiring rates pi and P2- Dark 
grey regions correspond to fragmentation with respect to A. 
The left panel shows the case pz = (pz = 0.5). This corre- 
sponds to an uncoupled system: the critical rewiring rates for 
pi and P2 are the same as for the two-state voter model, p c . 
In the right panel, pz is maximal (p3 = 0). Here, the active 
link density in the active cluster leads to an extension of the 
fragmentation region. Black lines represent analytical results. 
N — 10000, {k)—4, averaged over 20 realizations. 
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rates. As there is no analytical expression for p{jp) in 
the whole p-range ([26 ) the long-term behavior can only 
be predicted with certainty in regions I and III of the 
fragmentation diagram. 

In summary, evaluating the partial fragmentation con- 
dition with respect to A, i.e. A(J) = 0, leads to a phase 
diagram as shown in Fig. [3j where three different regions 
can be distinguished: In regions I and III partial frag- 
mentation occurs or is avoided regardless of P3, whereas 
in region II partial fragmentation depends on p$ and con- 
sequently on the setting of the related rewiring rate p$. 
We found that these results are in very good agreement 
with data obtained from agent-based simulation of large 
networks (Fig. 

Until now, we studied partial fragmentation with re- 
spect to one specific state (state A). In order to predict 
the final state of the whole system, one has to analyze the 
corresponding partial fragmentation diagrams for each of 
the three states. Let us assume p\ > p 2 > P3- Then, 
there are four cases to distinguish (see Fig.[5|: 

1- P3 > Pc 

2. p 3 < p c and p 2 > p c 

3. P2 < Pc and p 1 > p c 

4. pi < p c . 

In case 1) full fragmentation is reached, because all points 
P = (p 2 ,pi), Q = (P3,pi) and R = (p 2 ,P3) lie in the 
region III of their respective diagrams. In case 2) P lies 
in III and Q and R in I. Thus, partial fragmentation with 
respect to A occurs, while B- and C-nodes form an active 
cluster and in a finite system eventually reach consensus. 



FIG. 5. Fragmentation diagrams for partial fragmentation 
with respect to state A, C or B for p\ > P2 > Pa- There 
are four different cases, according to the conditions given on 
the left-hand side of the chart. The positions of the points 
P = {p2,Pi),Q = (p3,pi),R — (f>2,f>3) indicate for each case 
whether partial fragmentation is reached for the respective 
state. The final state of the whole system can be deduced 
from the outcomes for all three states and is given on the 
right-hand side of the chart. For the first, second and forth 
case the final state can be predicted without ambiguity. In 
the third case, either partial fragmentation with respect to 
A or consensus can be reached, depending on the specific 
parameter setting (the two possibilities are indicated by two 
different symbols o and x in the first diagram and a solid 
and a dashed line in the second diagram of the third row). It 
can be seen that partial fragmentation is only possible with 
respect to the state with the largest rewiring rates (here A). 



In case 3) the point P lies cither in region I or II of the 
corresponding fragmentation diagram for A, whereas Q 
is always in region I, because p 2 > Pz and P3 > p 2 - This 
means that in this case either partial fragmentation with 
respect to A or consensus is reached depending on the 
specific values of pi,p 2 ,p 3 . In case 4) P, Q and R lie 
in the region I of their respective diagrams and global 
consensus is reached. Note that this shows that partial 
fragmentation can only occur with respect to that state, 
which is connected via the largest rewiring rates to the 
two other states in the state network. 

In summary, we showed that in the three-state voter 
model either consensus, partial fragmentation or full frag- 
mentation occurs. Full fragmentation is only reached 
when all rewiring rates exceed p c . Analyzing the phase 
diagram with respect to the state which is connected by 
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the largest rewiring rates to the other states suffices for 
the prediction of the final state of the whole system. For 
quantitative predictions in region II of the diagram the 
active link density corresponding to the lowest rewiring 
rate has to be known. Qualitatively, one can say that 
if partial fragmentation occurs, then with respect to the 
"most different" state. 



IV. FRAGMENTATION TRANSITIONS IN A 
G-STATE VOTER MODEL 

Let us now consider a general system of G states. In 
contrast to the previous system partial fragmentations 
can also occur with respect to a group of states. A general 
multi-state network can thus fragment into several active 
components. Let us therefore calculate the condition for 
a system to fragment into two components containing s 
and G — s states, respectively (see Fig. [6]). This is in prin- 
ciple no restriction, as a fragmentation into more than 
two components can be treated as a fragmentation into 
two components where the active components in their 
turn fragment. 

For clarity we only use one level of indices from now 
on: we write Djj and instead of D Pij . and X Pi:j . 
Furthermore, in order to distinguish indices which re- 
fer to one component from those refering to the other 
component we use indices i G {l,...,s} for the com- 
ponent with s states and indices i 6 {1, . . . , s} for the 
component with s — G — s states. For example, the 
inter-component rewiring rates are then denoted as pa 
and intra-component active link densities as p^ and p^, 
respectively. 

In analogy to the treatment of the three-state model, 
we consider a situation where the two clusters are almost 
fragmented. We then determine the evolution equations 
for a set of active motifs connecting the two components. 
In the three-state case these were of two types, AB- and 
AC-fans, which led to a Jacobian of 2 x 2 matrix-valued 
entries and a fragmentation condition which was a func- 
tion of the rewiring rates p\ and P2 and the active link 
density p^. In the general case the Jacobian contains 
ss x ss matrix-valued entries, according to the number 
of inter-component links in the state-network (see Fig. [6| 
and the fragmentation condition is a function of all inter- 
component rewiring rates {pa} and all intra-component 
active link densities {pij} and {ptj}- 

Following the same procedure as for the three-state 
system, one finds that the general Jacobian exhibits a 
block-structure of s x s submatrices, 



J {Pii j Pij ) Pij ) — 

(£^(pij,pij) £i(pi 2 ) 

6(p2l) A 2 (p»j,P2j) 

V ?»(p«i) 



(5) 



£o-i) (P(«-i)») 

6(Ps( S -l)) &s{piJ,Psj) ) 



where Aj and are matrices of s x s matrix-valued en- 
tries, 



X «(p2l) t)i2(p2i,Pi]) 



\ Xi £ (psl) 

and 

6 (pij) 



Xjl(pis) \ 

Xfa-l( p(s-l)s ) 
X ia -l(p s ( a _ 1 )) T>is{psj_, pij) ) 



/X<l(p i; 



\ 



X-2(pij) 







V o ••• o xisipij)/ 

Here, we introduced the abbreviation 

S s 

i=hi¥i 3=1) 

The matrices D and X for k = 3 were already given in 
([2]) and (J3J) and the matrix X' for k = 3 is, 



'0 i(l - Pii ) Pi 



X-'u(Pij) =00 







v o Ki-ftiK-. 



(6) 



The latter matrix appears for fragmentations where both 
of the fragmenting components are active, i.e. for 1 < s < 
G-l. 

Note that every matrix-valued row of the Jacobian 
sums to Djj and refers to one specific type of inter- 
component link with rewiring rate pa in the state- 
network. One such row thus represents the transitions for 
one motif-type. Entries D on the diagonal capture the 
creation of motifs of the same type, while off-diagonal en- 
tries X and X' denote transitions to different motif-types, 
which arise from the intra-component link densities p^ 
and p^, respectively. For example, the entries in one 
row describing the transitions of (7152-fans depend on the 
rewiring rate between the states g± and g2, which is P12, 
the active link densities between g\ and all other states 
in the first component, pi x , with x G {2, . . . , s}; and the 
active link densities between 32 and all remaining states 
in the second component, p2 y , with y G {1,3, . . . ,s}. 

In analogy to the three-state model the active link den- 
sities pij entering in the Jacobian relate to the global 

active link density p£p as 



Pij 



[9i9j] 
k[9i] 



G [gigj] 
2 L 



G (G) 

2 h ■ 



(7) 



This holds analogously for pij. 

Stability analysis of the general Jacobian in ^ is in 
principle possible, but leads to a fragmentation condition 
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ces Mjj, such that 



FIG. 6. schematic representation of a state-network with G 
states. Ellipses illustrate a fragmentation into two (possibly) 
active clusters of s and G — s states. Dashed lines correspond 
to inter-cluster links with rewiring rates {pu}, connecting ev- 
ery state in one cluster with every state in the other cluster. 
The number of inter-cluster links, s(G — s), determines the 
dimension of the (matrix-valued) Jacobian. Solid lines corre- 
spond to intra-cluster links within both clusters with active 
link densities {pij} and {pij}, respectively. 



which depends directly or indirectly (through the active 
link densities) on all G{G — l)/2 different rewiring rates. 
In contrast to the estimation of the active link density in 
a two-state system, in a multi-state system the active link 
density of a certain link-type does not only depend on the 
rewiring rate of that specific link-type, but also on the 
rewiring rates and active link densities of the neighbor- 
ing links in the state-network. Inferring the link densities 
analytically from the rewiring rates is presently an un- 
solved challenge. So, even for given rewiring rates it is 
in general not possible to make quantitative predictions 
about fragmentation thresholds. Nevertheless, the struc- 
ture of the Jacobian allows for qualitative predictions, 
which will be shown in the next section. 
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as generalized, matrix-valued rowsums of M. Then, the 
following inequality holds [28] : 



A(minS fe ) < A(M) < A(maxS fe ) 

k k 



(8) 



The expressions min^ and max^ have to be understood 
element-wise, i.e. the matrix min^- is the matrix which 
is obtained when we take element-wise the minimum over 
all Sfc and analogously for the maximum. 

In the following we apply the theorem quoted above to 
the Jacobian given in ([5| . This is possible because J can 
be written as J = T 1, where T is a nonnegative irre- 
ducible matrix and 1 is the identity matrix of appropriate 
dimension. We will consider two different partitions. 

First, let us consider a partition (PI) of the Jacobian 
into ss submatrices. Then, the matrix-valued rowsums 
corresponding to this partition yield 

1 s 

= p^) + Xii(pij) + x'ii(pij) 



k = 1, . . . , ss. 



The matrices Djj only depend on pa and it can be seen 
from Q that all non-constant entries in D D increase with 
decreasing p. Therefore, we get for the upper and lower 
bounds of A(J) 



A(D Pmax )<A(J)<A(D Pmm ), 



(9) 



where 



V. REDUCTION PRINCIPLES FOR SPECIAL 
STATE-NETWORK TOPOLOGIES 

In this section we use theorems about upper and lower 
bounds of the largest eigenvalue A(M) of a nonnega- 
tive irreducible matrix M. The well-known Frobenius 
inequality states 

minSfe < A(M) < max.Sk, 

k k 

where St is the rowsum of the i-th row of M. A general- 
ization of the above inequality for a partitioned nonneg- 
ative irreducible square matrix M is given in |28j . Let us 
assume that M can be partitioned into square submatri- 



Pmax = max pu, p min = minpu. 

From ([9]) we deduce the following statements: 

• When all inter-cluster rewiring rates pu are below 
the threshold p c , no fragmentation occurs, because 
A(D Pmax ) > 0. 

• When all inter-cluster rewiring rates pu exceed 
the threshold p c , fragmentation occurs, because 
A(D Pmin ) < 0. 

• If Pmin = Pmax, necessarily all inter-cluster rewiring 
rates must be equal. In that case, the fragmenta- 
tion condition is the classical condition of the two- 
state voter model, A(J) = A(D p ) = 0, which yields 
the critical rewiring rate p c . 
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FIG. 7. Numerical test of the reduction principle. We run 
simulations for a 5-state system assigning random values to 
the rewiring rates r; (inset). Plotted is the density of inter- 
cluster links (dashed links) for three different values of p. 
For each p-value we run ten independent simulations with 
randomly chosen ri and compare the corresponding inter- 
cluster link density to the active link density of a two-state 
model with the same rewiring rate and a ratio 2:3 for the 
number of nodes in opposite states. It can be seen that 
the inter-cluster density does not depend on n and that the 
steady-state value equals the corresponding inter-cluster den- 
sity (active link density) in the original two-state voter model. 
N = 10000, (k) = 4. 



The first two results represent an intuitive generalization 
of our findings for the three-state case. The last result 
implies that if all inter-cluster rewiring rates are equal 
then the value of these rewiring rates, p, is the only pa- 
rameter parameter on which the fragmentation transition 
depends. In this case a precise analytical estimation of 
the fragmentation point is possible because the active 
link densities arising from the intra-cluster links, do not 
enter. Furthermore note that this result is independent of 
the number of opinions. This implies that in the special 
case of equal inter-cluster rewiring rates systems of any 
size behave identically to a properly-initialized adaptive 
two-state voter model. 

We test the latter result in a five-state system, con- 
sidering a fragmentation into two components of 2 and 
3 states, respectively (inset in Fig. [7]). Simulations show 
that for randomly chosen rewiring rates within the two 
clusters, the inter-cluster link density reaches the same 
steady-state value (Fig. [7]). A further comparison shows, 
that the behavior of a five-state model closely matches 
that of a two-state model. 



Now we consider another partition, (P2), of the Jaco- 
bian, which is a partition into s submatrices. Then, the 
corresponding generalized rowsums yield 



S T 2) = A, (//,,./,,,) + &(f*j) 

/Dii(py) Xii(pi2) 
Xi2(p2i) Di2(p2j) 



Xjl(pls) \ 



Xi s _l(p( s _ 1 ) s 
Xi a-l (p s ( a _ 1 )) Dis(psj) / 



where 



~Dii(pii) = D.i - Xij(py). 

First, we observe that every matrix corresponds 
to a partial fragmentation in a system of s + 1 states. 
More precisely, the set {S- ; } describes a collection of 
s single-state-fragmentations where for every i a single 
state gi is taken separately from the s-cluster. This single 
state (i.e. now s = 1) then forms the first component of 
the partial fragmentation, while the second component is 
given by the whole s-cluster. 

Now, building the element-wise extrema of {S,[ P2 ''} 
means to compare all the single-state-fragmentations by 
comparing every matrix-entry of the corresponding gen- 
eralized rowsums. Taking mini S^ P2 ^ (max.; s[ F2 ^ ) yields 
therefore in evere matrix-entry the minimum (maximum) 
value, i.e. that one which comprises the maximal (min- 
imal) rewiring rate. The resulting matrix corresponds 
to a partial fragmentation with respect to a single state 
where the inter-cluster rewiring rates are chosen extremal 
according to the described comparison. We will refer to 
such a system as bounding system (see Fig. [8] for exem- 
plary bounding systems). 

For a partition of type (P2) the leading eigenvalue of 
the general Jacobian satisfies 



A(minSf 2) ) < A (J) < A(maxSf 2) 



(10) 



The lower bound corresponds to the fragmentation of a 
system where the largest inter-cluster rewiring rate of 
each state in the second component is connected to a 
single state. The upper bound corresponds to the frag- 
mentation of a system where the smallest inter-cluster 
rewiring rate of each state in the second component is 
connected to a single state. 

As in the Jacobian in ([5| the matrices X and X' can 
be interchanged, one can consider a corresponding parti- 
tion where the second component is reduced to a single 
state, i.e. s — 1 and the first component remains as a 
whole. This leads to a different set of bounding systems, 
as illustrated in Fig. [8] a) and b). 

From ( 10 ) we can draw the following conclusions: 
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FIG. 8. Schematic representation of the generalized Frobe- 
nius inequalities using the example of a five-state system and 
a specific partial fragmentation into a 2-state- and a 3-state- 
cluster. The inequalities in a) and b) represent two different 
sets of bounding systems obtained from a partition of type 
(P2), as described in the text. In a) the bounding systems 
are three-state systems, obtained by reducing the 3-cluster to 
a single state. The remaining inter-cluster rewiring rates are 
the respective maximum (minimum) values of the rewiring 
rates between each state in the 2-cluster and all states in 
the 3-cluster. In b) the bounding systems are four-state sys- 
tems, obtained by reducing the 2-cluster to a single state. 
The remaining inter-cluster rewiring rates are the maximum 
(minimum) values of the rewiring rates between each state 
in the 3-cluster and all states in the 2-cluster. In c) a special 
case is shown where the upper and lower three-state bounding 
systems coincide. In this case, the leading eigenvalue of the 
five-state system equals the leading eigenvalue of the associ- 
ated three-state system and the considered fragmentation of 
the original system is fully captured by the lower-dimensional 
one. 

• If the lower bounding system does not fragment 
(A(mirij S,- P2 ^) > 0) the original system does not 
fragment. 

• If the upper bounding system fragments 
(A(maxi S- P2 ^) < 0) the original system frag- 
ments. 

• If upper and lower bounding systems are the same 
(mini s' P2 ^ = maxi sj P2 ^) the fragmentation of the 
original system is exactly captured by the bound- 
ing system, i.e. the full Jacobian J reduces to the 
Jacobian of the bounding system. 



The latter case is realized if every state in one compo- 
nent is connected via equal rewiring rates to every state 
in the other component (see Fig.[8]c)). For state-network 
topologies which display this property the dimension of 
the Jacobian reduces significantly and thus the fragmen- 
tation condition becomes much more tractable. 

In summary, the results from the second partition show 
that for the leading eigenvalue of a Jacobian, correspond- 
ing to a partial fragmentation into two active clusters, 
upper and lower bounds can be given, which correspond 
to single-state-fragmentations in (properly constructed) 
lower-dimensional systems. In particular, the leading 
eigenvalue of the full Jacobian can be exactly calculated 
as the leading eigenvalue of a lower-dimensional Jacobian 
if special state-network topologies are given. Otherwise, 
when such a reduction is not possible, the bounding sys- 
tems provide necessary conditions for a partial fragmen- 
tation to occur. So, calculating fragmentation conditions 
for the much simpler bounding systems in some cases 
suffices to predict the ocurrence or absence of fragmen- 
tations in the full system. 



VI. CONCLUSION 

In the present paper we extended recent work on the 
adaptive two-state voter models to a family of multi-state 
models. For the three-state model our analysis revealed 
a phase diagram in which three distinct types of long- 
term behavior are observed. Depending on the param- 
eters the system either approaches a consensus state, a 
partially fragmented state ultimately leading to two sur- 
viving opinions or a fully fragmented state in which all 
three opinions survive. 

In a general scenario with an arbitrary number of states 
making precise predictions is more difficult. In particu- 
lar, the computation of fragmentation points generally 
requires the estimation of active link densities inside the 
clusters between which the fragmentation occurs. By ex- 
ploiting the specific structure of transition rates in the 
system, one can nevertheless gain analytical insights into 
the fragmentation dynamics. For example we identified a 
class of special cases in which adaptive multi-state voter 
models exactly recover the behavior of the adaptive two- 
state voter model. 

While the ultimate goal of understanding opinion for- 
mation in the human population is still far away, the 
present progress shows that analytical understanding can 
be pushed to more complex models. Recent studies have 
shown that already today variants of the voter model can 
be tested in experiments with swarming animals. An im- 
portant goal for the future is to continue the refinement of 
models and analysis techniques in order to describe real- 
world situations. We hope that the approach presented 
here will make a contribution to this ongoing process. 
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VII. APPENDIX 

We use a moment closure approach [29] for the calcula- 
tion of the maximal active link density in an active cluster 
of s states in a G-state system. The evolution equation 
for the number of active links of type xy is given by 



the number of nodes in the active component. The total 
number of links in the active component I is then given 
by 

I = ST] + 



[xy] = -[xy] + ^(^2[xzy]+^2[xzy}-Y^[xyz}-^2[zxy}), and © yiclds 



z^x z^y z^y z^x 

assuming p,j = for all rewiring rates within the active 
cluster. Then, using the pair-approximation, we get for 
the steady-state, 

r i = 1 f 2 HM + 2 NHyy] + 2 HN 
9 ^ 1 ' [y] N 



2 V ^ ~ [z 



N] [yz] _ [zx] [xy] \ 
^ [y] ^ [x] )' 

z^y lal z^x L 1 



(11) 



C = ~s(s - 2) 



n 

'■2 



Cv C 2 

n n 

CI C 2 
-s(s-2) + 2— -2— s(s-l). 



Using kn = 21, we get for the maximum active link den- 
sity in a cluster of s states: 

_ s( a -l)C _ (s-l)(fc-l) 
Pmax ~ 21 ~ sk 



For equal distribution of states we can write [a;] = 
n/sVx, [xy] = (\/(x,y) and [xx] — r/\fx, where n denotes 



which gives p max = (k — l)/(2k) for s = 2, as provided 
in the text. 
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